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Abstract 

Factor H-Binding Protein (fHbp) is an outer membrane protein antigen included in two novel meningococcal group B 
vaccines and, as sucli, is an important typing target. Approximately 50% of meningococcal disease cases in England and 
Wales are confirmed using real-time PCR on non-culture clinical specimens only. Protocols for typing fHbp from this subset 
of cases have not yet been established. Here we present a nested PCR-based assay designed to amplify and sequence fHbp 
from non-culture clinical specimens. From analytical sensitivity experiments carried out using diluted DNA extracts, an 
estimated analytical sensitivity limit of 6 fg/\xL of DNA (<3 genome copies/nL) was calculated. The sensitivity of the assay 
was shown to be comparable to the ctr/A-directed real-time PCR assay currently used to confirm invasive disease diagnoses 
from submitted clinical specimens. A panel of 96 diverse, patient-matched clinical specimen/isolate pairs from invasive 
disease cases was used to illustrate the breadth of strain coverage for the assay. All fHbp alleles sequenced from the isolates 
matched those derived from previous whole genome analyses. The first-round PCR primer binding sites are highly 
conserved, however an exceptional second-round PCR primer site mismatch in one validation isolate prevented 
amplification. In this case, amplification from the corresponding clinical specimen was achieved, suggesting that the use of 
a nested PCR procedure may compensate for any minor mismatches in round-two primer sites. The assay was successful at 
typing 91/96 (94.8%) of the non-culture clinical specimens in this study and exhibits sufficient sensitivity to type fHbp from 
the vast majority of non-culture clinical specimens received by the Meningococcal Reference Unit, Public Health England. 

Citation: Clark SA, Lucidarme J, Newbold LS, Borrow R (2014) Genotypic Analysis of Meningococcal Factor H-Binding Protein from Non-Culture Clinical 
Specimens. PLoS ONE 9(2): e89921. doi:10.1371/journal.pone.0089921 

Editor: Baochuan Lin, Naval Research Laboratory, United States of America 

Received November 27, 2013; Accepted January 29, 2014; Publislied February 24, 2014 

Copyright: © 2014 Clark et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was funded by Pfizer Vaccine Research, 401 N. Middletown Rd., Pearl River, New York, United States. The funders had no role in study design, 
data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: This work was funded by Pfizer Vaccine Research, 401 N. Middletown Rd., Pearl River, New York, United States. RB is a member of the 
PLOS ONE editorial board. SAC and LSN have performed contract research on behalf of Public Health England for Pfizer Vaccine Research. JL and LSN have 
performed contract research on behalf of Public Health England for Novartis Vaccines and Diagnostics. RB has performed contract research on behalf of Public 
Health England for Baxter Biosciences, GSK, Merck, Novartis Vaccines and Diagnostics, Pfizer Vaccine Research, Sanofi Pasteur and Sanofi Pasteur MSD. This does 
not alter the authors' adherence to PLOS ONE policies on sharing data and materials. 

* E-mail: stephen.clark@phe.gov.uk 



Introduction 

Invasive meningococcal disease (IMD) represents a significant 
threat to the public health of communities globally. Immunisation 
against the causative organism, Neisseria meningitidis, remains the 
most effective means of curtailing the prevalence of this deadly 
disease. Vaccines containing peptide-conjugated capsular polysac- 
charides have been successful in reducing disease caused by 
capsular groups A, C, W and Y [1-4]. Disease caused by group B 
meningococci, however, remains an ever present problem and the 
development of vaccines against this subset has proven to be much 
more challenging. Structural similarities between group B 
polysaccharide and human foetal antigens inhibit the immunoge- 
nicity of polysaccharide-based group-B vaccines [5]. Furthermore, 
formulations based upon group B outer membrane vesicles 
(OMVs) are yet to offer expansive immune protection due to 
immunodominance of the heterogeneous Porin A (PorA) antigen 
[6] . In a change of approach, two group B vaccines containing 
subcapsular, outer membrane protein antigens have recendy been 
developed and currently represent the most promising strategy for 
preventing diverse group B disease [7] . Factor H-Binding Protein 
(fHbp) is an antigen included in both of these vaccines. 



FHbp (previously known as LP2086 and GNA1870) is a surface- 
exposed lipoprotein composed of two domains (a P-sheet and a |3- 
barrel) with an N-terminal flexible chain that anchors it to the 
outer membrane [8,9]. The antigen was identified over a decade 
ago by two separate groups and is perhaps one of the most studied 
subcapsular vaccine candidates [10,11]. The antigen is found 
among virtually all invasive isolates and enhances survival in vivo 
through its ability to bind human complement factor H, a down- 
regulator of the host alternative complement pathway [12]. 

At the time of writing, alleles coding for over 600 variants of the 
fHbp peptide have been sequenced and deposited into the 
dedicated PubMLST database (http://pubmlst.org/neisseria/ 
fHbp) [13]. These peptide variants can be categorised using two 
distinct nomenclatures based upon amino acid sequence identity. 
One system divides variants into two subfamilies (A and B) [10] 
and the other into three variant groups (1, 2 and 3) [1 1]. Subfamily 
B (sB) corresponds to variant group 1 (vl) whilst subfamily A (sA) 
encompasses variant groups 2 and 3 (v2 and v3). The PubMLST 
database (hosted by the University of Oxford, UK) seeks to 
reconcile these differences by assigning arbitrary, sequential 
numbers to alleles and peptides as they are submitted. These are 
then annotated with the corresponding designations from the two 
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systems described above [14]. The PubMLST nomenclature will 
be primarily referenced within this article. 

Whilst considerable amino acid sequence variation occurs in 
fHbp across invasive populations, several studies have illustrated 
the cross protective nature of the antigen whereby one fHbp 
variant can elicit immune protection against strains bearing 
variants of the same subfamily [8, 1 0, 1 1 , 1 5] . This attribute of 
subfamily-specific cross protection is exploited in an investigational 
recombinant vaccine (rLP2086), developed by Pfizer. The vaccine 
contains two lipidated fHbp peptide variants, peptide 45 (sA/v3) 
and peptide 55 (sB/vl), and trials in adults have shown it to elicit 
serum bactericidal antibodies (SBA) against diverse group B 
isolates [16,17]. In a parallel effort, a multi-component group B 
vaccine, also containing fHbp, has recently been licenced for use in 
children and adults in the European Union and Australia [18,19]. 
The formulation contains one fHbp variant, peptide 1 (sB/vl), as 
part of a fusion protein among other antigens [20]. 

The potential use of these vaccines in the UK and elsewhere 
emphasises the need for comprehensive enhanced surveillance of 
disease-causing organisms. These activities aid vaccine coverage 
predictions and wiU serve to monitor the ongoing impact of the 
vaccines should they be utilised. In England and Wales, Pubhc 
Health England's Meningococcal Reference Unit (PHE, MRU) 
provides species confirmation and routine characterisation of IMD 
isolates. Since July 2010, whole genome sequence analysis has 
been adopted for all invasive isolates received by the MRU. To 
date, the genomes of over 900 UK isolates (epidemiological years 
2010/201 1 and 201 1/2012) have been sequenced, annotated and 
made freely available in the Meningitis Research Foundation's 
online Meningococcus Genome Library (MGL, http://www. 
meningitis.org/current-projects/genome). In recent years, howev- 
er, approximately 50% of IMD cases in England and Wales were 
confirmed solely by real-time PGR testing of clinical specimens 
from which no isolate could be obtained [21]. Limiting factors 
such as low DNA concentration and interference from human 
DNA have prevented the adoption of whole genome analysis from 
these 'non-culture' specimens. 

PCR-based assays designed for MultUocus Sequence Typing 
(MLST) and typing of PorA loci from clinical specimens have been 
developed [22,23]. An assay to characterise fHbp in this way is, 
however, not currentiy available. Gonsequentiy, the variation of 
fHbp among this large proportion of clinical cases remains 
unknown. This represents a significant epidemiological knowledge 
gap and may influence coverage estimations of fHbp-containing 
vaccines. In order to complete the picture of fHbp distribution in 
England and Wales, a PCR-based assay has been developed to 
sequence fHbp from non-culture clinical specimens containing 
potentially low-copy number meningococcal DNA. This work 
describes the assay and the activities carried out to illustrate its 
sensitivity to the target template and its coverage of diverse 
meningococcal strains. 

Methods 

Ethics Statement 

As this study involved the use of anonymised specimens that 
were surplus to cUnical need, ethical approval was not required in 
accordance with National Research Ethics Service (NRES) 
guidelines. 

Isolates and Clinical Specimens 

Analytical sensitivity. Seven invasive meningococcal iso- 
lates, representing the predominant clonal complexes and fHbp 



variant groups, were selected to assess the analytical sensitivity of 
the assay (Table SI). 

Validation panel. A validation panel comprising 96 clinical 
isolates was compiled from disease cases confirmed by both 
organism isolation and PGR on a clinical sample (2010-2012). All 
isolates selected ha\'e previously undergone whole genome analysis 
and are included in the MGL. Each isolate was paired with its 
corresponding patient-matched clinical specimen. Selection of the 
isolate/specimen pairs was based upon MLST and fHbp data of 
the cultured isolate. The panel was selected to represent diverse 
clonal complexes (GG), serjuence t)'pes (ST) and fHbp variants, 
approximating the distribution of GGs within England and Wales 
(Table SI). The diversity of these specimens was confirmed by 
mapping the respective STs against their corresponding GGs as 
determined by eBURST V3 analysis [24] of the entire MLST 
dataset (PubMLST database, 30/12/2013, Fig. SI). To improve 
tiie resolution of the ST- 11 GC, ribosomal MLST (rMLST) and 
SplitsTree4 analysis (www.splitstree.org) was performed on avail- 
able ST-1 1 GG genomes, against which the selected panel isolates 
were mapped [25]. 

Neisseria lactamica isolates. In addition to the vahdation 
panel, a selection of six diverse N. lactamica isolates was made based 
upon the distribution of the STs within the known population of 
the organism [26]. Isolates from distinct STs and GGs were chosen 
(Table SI). Primer complementarity was assc'ss(;d using an 
alignment of sequences from 16 jV. lactamka isolates submitted to 
die PubMLST database (Table S2). 

DNA Extraction 

All DNA extractions were performed using the Qiagen DNeasy 
Blood and Tissue Kit. 

DNA extraction from clinical specimens was carried out in 
accordance with Qiagen DNeasy Blood and Tissue Handbook 
(July 2006). A volume of 100 |iL of each clinical specimen was 
used where possible. If less than 100 |iL was available, the 
maximum volume available was used and the remainder was 
made up using phosphate-buffered saline (PBS). Extracts from 
clinical specimens were eluted twice using 50 nL buffer AE per 
elution. 

DNA extraction from cultured isolates was adapted from the 
manufacturer's gram-negative bacteria protocol (Qiagen Blood 
and Tissue Handbook, July 2006). All isolate extracts were eluted 
in 2x75 H-L buffer AE. Neisseria lactamica isolate extracts prepared 
previously were used [26]. All DNA extracts were stored at 4°G. 

Analytical Sensitivity 

Extracted DNA from each isolate was quantified using Qubit 
2.0 fluorometer (Life Technologies, USA) following the manufac- 
turers' protocol and adjusted to 6 ng/^iL. A series of eight 10-fold 
dilutions was carried out to an estimated DNA concentration of 
60 ag/ixL (Table 1). 

Semi-quantitative cycle threshold (Gt) values were obtained 
using the current in-house ctrA-dircctcd TaqMan assay, which is 
an adaptation of a previously published protocol [27]. The assay 
was performed using the ABI 7500 Fast Real-Time PGR System 
(Life Technologies, USA) and features cyanine dye (Gy5) probes 
targeting sequences within the ctrA capsule syntli('sis gene. 
TaqMan assay plates with pre-lyophUised primer/probes (Life 
Technologies, USA) were reconstituted in molecular grade water 
before 5 (i,L of the diluted DNA extracts was added. The Gt value 
is given as the cycle-number at which the fluorescent signal meets 
a calculated background threshold, confirming the presence and 
amplification of the target sequence. A Gt value of >45 is 
currentiy reported as negative. 
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Table 1. Serial dilutions of meningococcal extracts, corresponding estimated DNA concentration and calculated genome copy 
number (estimate based upon a 2.2 IVlbp genome). 



Dilution 


DNA Concentration (per \iL) 


Calculated Single Genome Copies (per ^L) 


Starting Cone. 


6 ng 


2.53x10'' 


10-' 


600 pg 


2.53x10= 


10-^ 


60 pg 


2.53x10" 


10-^ 


6 pg 


2530 


10-" 


600 fg 


253 


10-= 


60 fg 


25.3 


10"^ 


6 fg 


2.53 


10-^ 


600 ag 


0.253 


IQ-" 


60 ag 


0.0253 


doi:l 0.1 371 /journal.pone.0089921 .tool 



Finalised Assay Protocol 

Detailed descriptions of the assay development (including 
primer design and optimisation) are included as Method SI. 

Nested PGR. To attain the sensitivity and specificity required 
for amplification from non-culture specimens, a nested-PCR 
protocol, utilising two consecutive PGR rounds, was adopted. 
PGR primers are listed in Table 2. All PGR reagents and volumes 
are listed in Table 3. Round one PGR reactions were performed in 
a ,50 |J,L reaction containing 10 [iL of extracted template DNA. 
Round two reactions were performed using a 25 |iL reaction 
containing 2 |J.L of round one product. Negative controls were 
prepared by replacing DNA extract with 10 |J,L of molecular- 
grade water in the round one reaction. This was then transferred 
to PGR round two along with experimental products. 

Both PGR rounds use the same temperature and thermocycling 
conditions /times. An initial activation step of 95°G for 15 minutes 
was followed by 35 (round one) or 30 (round two) cycles of 
denaturation at 95°G for 30 seconds, annealing at 63°G for 30 
seconds and extension at 72°G for 80 seconds. These were 
followed by a final extension step of 72°G for 7 minutes. Following 
PGR round two, amplification was confirmed by 2% agarose gel 
electrophoresis. 

Successfully amplified products were purified prior to sequence 
analysis using ExoSap IT PGR product cleanup kit (Affymetrix, 
USA). The purified product was diluted in molecular grade water 
based upon visual assessment of gel band intensity. A 1 in 4 
dilution was performed for concentrated products exhibiting 
relatively strong gel bands, whilst products with average intensity 



gel bands were diluted 
diluted 1 in 2. 



in 3. Products with weak gel bands were 



Sanger Sequencing 

Sequence analyses were performed using Big Dye Terminator 
v3. 1 cycle sequencing kit (Life Technologies, USA). Big Dye 
terminator reaction mix was used in a 1/8''' reaction according to 
the manufacturer's protocol (2002). Products were cleaned by 
sodium acetate/ethanol precipitation and resuspended in 15 (il of 
Hi-Di Formamide (Life Technologies, USA) before analysis using 
an ABI 3130x1 Genetic Analyser. Sequence assembly and analysis 
was carried out using Sequencher v4.7 (Gene Godes Gorporation, 
Ann Arbor, MI USA). 

Results 

Analytical Sensitivity 

To assess the sensitivity of the PGR protocol to the target 
template, DNA was extracted from seven diverse meningococcal 
isolates and diluted in series to an approximate concentration of 60 
ag/|a,L (0.025 genome copies (GG)/|xL) (Table 1). A volume of 
5 |J.L of extract dilution was initially used in a 25 |J.L round one 
PGR reaction. PGR products of the appropriate size were visible 
on the agarose gel after the first PGR round for all extracts down 
to a DNA concentration of 60 fg/|iL (25.3 GG/|.iL). For PGR 
round two, varying amounts of round one product (1 (iL, 2 |J.L 
and 5 |J.L) were added to a 25 |iL final reaction volume. After the 
second PGR round, detectable amplification occurred in all 



Table 2. PCR and 


sequencing primers. 








Primer ID 


Primer Use (direction) 


Sequence (5' to 3') 


Reference 


1869-2F 


Round One PCR (Fwd) 


GAAGAAATCGTCGAAGGCATCAAAC 


126,29] 


1871 Rait 


Round One PCR (Rev) 


ATGCCGATACGCAGTCC(G/C)GTAAAC 


[26,29] 


fHbpRd2F 


Round Two PCR (Fwd) 


GTTATGCCAAGGGCGAATTGAACC 


n/a 


fHbpRdlR 


Round Two PCR (Rev) 


GTGCGGATTTCCGGCAG(G/A)ATCA 


n/a 


gnal870F 


Sequencing (Fwd) 


TGACCTGCCTCATTGATGC 


[30] 


fHbpseqRI 


Sequencing (Rev) 


AGGACGGG(G/A)CGGTT(G/A)AAATC 


n/a 


gnal870S2 


Sequencing (Fwd) 


CAAATCGAAGTGGACGGGCAG 


[30] 


gnal870S3 


Sequencing (Rev) 


TGTTCGATTTTGCCGTTTCCCTG 


[30] 


doi:l 0.1 371 /journal.pone.0089921 .t002 
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Table 3. Standard nested-PCR mastermix reagent volumes. 







PCR Round 


Reagent 


Volume per Reaction (|iL) 


Round One 


lOX PCR Buffer 


5 




1869-2F (5 fiM) 


5 




1871Ralt (5 nM) 


5 


DNTP Mix 1 




HotStarTaq 


0.25 




Molecular-grade Water 


23.75 




DNA Extract 


10 


Round Two 


lOX PCR Buffer 


2.5 




fHbpRd2F (5 nM) 


2.5 




fHbpRd2R (5 (iM) 


2.5 




DNTP Mix 


0.5 




HotStarTaq 


0.125 




Molecular-grade Water 


14.875 




Round One Product 


2 



doi:l 0.1 371/journal.pone.0089921 .t003 



extracts to a DNA concentration of 6 fg/|xL (2.53 GC/|J.L). At 
600 ag/|j,L (<1 GC/|J,L), five of the seven extracts produced 
appropriate amplicons. No amplification was seen at 60 ag/ |a,L. 
The volume of round one product transferred into tlie round two 
reactions determined the intensity of tlie resultant bands, however, 
increasing the volume did not induce amplification in otherwise 
PCR-negative extracts. Using 2 |J.L of round one product was 
therefore deemed sufficient and this volume was used for all 
subsequent reactions. 

To determine the influence of DNA extract volume in the round 
one reaction, three extract volumes of the two highest dilutions 
(10"' and 10"") were added to round one reactions. The 
proportion of extract used in each reaction was kept constant: 

5 [xL/25 |xL, 10 nL/50 |aL and 20 jiL/lOO ^L. The amount of 
extract used in round one was shown be an important factor in the 
likelihood of gaining positive amplification. Figure 1 shows the 
round two products obtained following the respective round one 
reactions. The 600 ag/|j.L (10"') extracts from aU seven isolates 
produced appropriate amphcons when 20 |J,L of the diluted 
extract was used for the first round. Amplification occurred from 
only six and three of the seven extracts when 1 0 \xL and 5 |a,L of 
extract were used, respectively. PCR products were visible for four 
of the seven isolates when 20 |iL of the 60 ag/|j,L (10 ^) dilution 
was used. No amplification was seen in the other reactions at this 
dilution. 

In order to compare the DNA concentrations of these dilutions 
to routinely submitted specimens, the dilutions were applied to the 
in-house cto4-directed TaqMan real-time PCR assay. As expected, 
Ct values of the extracts steadily increased as the DNA 
concentration within each dilution decreased (Table 4). An 
average increase of approximately 4 cycles was observed for each 
successive extract dilution. Each extract was tested once up to 

6 fg/|J-L (10 *"). For dilutions of 10 ' and 10 the assay was 
repeated three times to counter any of the stochastic sampling 
effects observed previously. At 6 fg/ |J.L (the nested-PCR analytical 
sensitivity limit) an average Ct value of 37.86 was obtained. At 
10"', only one of the extracts produced a positive result and 
achieved a Ct value of 4 1 . No positive Ct value was obtained for 
any extract at the lowest dilution (10 "). 



Validation Panel 

A total of 96 diverse meningococcal isolates with corresponding 
non-culture clinical specimens were selected. DNA extracts from 
the clinical specimens were tested using the finalised nested-PCR 
protocol, however, the isolate extracts were tested using both PCR 
rounds independently. 

For PCR round one, all 96 isolates produced an appropriate 
amplicon. During PCR round two, aU but one isolate produced a 
PCR product (Table S3). Three attempts failed to produce an 
amplicon for isolate Mil 240189 for PCR round two. To 
determine primer site conservation, the JHbp allele sequence plus 
flanking regions for this isolate were aligned against the second- 
round PCR primers. Three base mismatches were found within 
the JHhpRd2R binding site (Figure 2), one of which was located at 
the 3' end from which extension occurs. Following a BLAST 
search of the 923 isolates within the MGL data, only one other 
isolate (Mil 240984) was found to feature these mismatches 
within this primer binding site (together representing 0.22% of 
MGL isolates). Sequencing of the 95 successfully amplified 
products revealed JHbp alleles matching the alleles assigned in 
the MGL from whole genomic data. Ml 1 240402 is SLnJHbp null 
isolate; no sequencing trace was obtained for this product. 

For the clinical specimens, 83 of the 96 (86.5%) produced an 
appropriate PCR product using the nested PCR procedure on the 
first attempt (Run 1, Table S3). PCR negative extracts were 
repeated in two subsequent runs yielding six further products 
(Runs 2 and 3, Table S3). The seven remaining PCR-negative 
extracts were tested at an increased extract/reaction volume 
(20 |J,L/100 |J.L) and PCR cycle numbers (round one increased by 
10 cycles to x45, round two increased by 5 cycles to x35. Run 4, 
Table S3). Five of these seven extracts stiU failed to amplify. These 
specimens were considered PCR negative and no further attempt 
at amplification was performed. After four attempts, therefore, 91 
of the 96 of the clinical specimens yielded PCR products (94.8%). 
Curiously, the specimen corresponding to isolate Mil 240189 
{tesLtunngJHbpRd2R site mismatch) amplified successfully using the 
nested PCR procedure. This may indicate that primer site 
mismatches can be overcome by using a preceding PCR round. 

All but two of the amplified products from specimens contained 
JHbp alleles matching those of their corresponding isolates. The 
alleles of the two incongruent specimens, Mil 922854 and Mil 
915746, both differed from their corresponding isolates' alleles by 
one base resulting in alleles not assigned in the PubMLST 
database. Re-amplification and sequencing of the specimens was 
performed to account for PCR copying errors that may have 
occurred. The resulting alleles matched the corresponding isolates 
and confirmed that errors were likely to have occurred early in the 
PCR cycle repeats. 

Neisseria Lactamica Isolates 

To provide insight into the conservation of the PCR primer sites 
beyond the species level, DNA extracts from 6 diverse jV. lactamica 
isolates [26] were tested in each PCR round individually (Table 
SI). Conservation of the primer site sequences in a related species 
could be indicative of functional constraint of the peptides for 
which they code and would increase confidence in the longevity of 
the primers. Recent work has confirmed that jV. lactamica do not 
possess JHbp, and that the allele for NLA18150 (a putative opacity 
protein) can be found in its place at this locus flanked by genes 
common to both species [26,28]. The PCR primers, targeting 
sequences within these flanking genes, should however produce an 
amplicon if the binding sites are conserved. For PCR round one, 
amplification occurred from all 6 jV. lactamica isolates. For PCR 



PLOS ONE I www.plosone.org 



4 



February 2014 | Volume 9 | Issue 2 | e89921 



Typing of fHbp from Non-Culture Clinical Specimens 



1 2 3 4 5 6 7 -ve 1 2 3 4 5 6 7 -ve ML 

5 |jL 

Extract in 

25 mL 
Reaction 

1 2 3 4 5 6 7 -ve 1 2 3 4 5 6 7 -ve ML 

10|jL 
Extract in 

50 [iL 
Reaction 

1 2 3 4 5 6 7 -ve 1 2 3 4 5 6 7 -ve ML 

20 mL 
Extract in 
100 \^L 
Reaction 



10-^ Extract Dilution 

Figure 1. Gel image of amplicons produced using varying extract volumes. Annotated gel image of round two PCR products from nested 
PCR using varying volumes of diluted DNA extracts in the preceding round one reaction. The products for all seven isolate extracts at the highest two 
dilutions are shown. Wells to the left of the image contain 10^' dilution extracts (600 ag/|iL) and wells adjacent to the molecular weight ladder 
contain 10"^ dilution extracts (60 ag/|iL). The extract/reaction volumes (|iL) are indicated on the right of the image. Extract numbers refer to the 
order in which the isolates are listed in Table SI. ML = Molecular ladder 
doi:1 0.1 371/journal.pone.0089921 .gOOl 

round two, however, amplification was only observed in 2 of the 6 
isolates. 

In order to determine the cause of the negative results, DNA 
sequences of NLA18150 and the flanking regions from 16 jV. 
lactamica genomes were downloaded from the PubMLST database 
and aligned with the PCR primer sequences. The genomic 
sequences were well conserved for all primer sites apart from 
JHhpRd2R. At this site, 4 of the 1 6 jV. lactamica sequences contained 



Table 4. Ct values generated using real-time TaqMan PCR assay. 





Isolate 


Ct Value 


per Extract Dilution 






















10^ 


10^ 


10^ 


10" 


10^ 


10" 


10 'a 


10 'b 


10 'c 


10 "a 


10"b 


10 'c 


M08 0240297 


16.8 


20.1 


25.0 


28.3 


32.1 


36.3 


41.0 


— ve 


— ve 


— ve 


— ve 


— ve 


M07 0240954 


17.4 


20.7 


25.3 


28.9 


32.3 


38.7 


— ve 


— ve 


— ve 


— ve 


— ve 


-ve 


M07 0241 036 


20.6 


25.8 


29.0 


33.0 


38.1 


41.0 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


M07 0241073 


17.0 


20.8 


25.4 


28.7 


32.9 


38.5 


— ve 


— ve 


— ve 


— ve 


— ve 


-ve 


M07 0240725 


17.4 


22.0 


25.2 


— ve 


32.3 


36.1 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


M08 0240113 


17.5 


22.1 


25.3 


27.6 


32.2 


36.5 


— ve 


— ve 


— ve 


— ve 


— ve 


-ve 


M08 0240032 


17.2 


21.8 


25.3 


28.4 


32.0 


37.9 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


Neg. Control 


-ve 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


— ve 


-ve 


Mean CT 


17.70 


21.85 


25.79 


29.15 


33.13 


37.86 


41.00 


n/a 


n/a 


n/a 


n/a 


n/a 



doi:1 0.1 371/journal.pone.0089921 .t004 
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10-s Extract Dilution 



the identical primer mismatches observed in the meningococcal 
validation isolate Mil 240189 (Figure 2). Five of the remaining 
sequences contained a single base substitution at position 9 of the 
primer (3 '-5'), however this is unlikely to affect amplification. 

To further assess the potential ability of the nested PCR design 
to compensate for round two primer site polymorphisms, 2 |J.L of 
round one product from all 6 .A^. lactamica isolates was added to 
25 |J.L round two reactions. The use of round one product as the 
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fHbpRd2R TGATYCTGCCGG2VAATCCGCAC 

(A) GAAC CCTA 

(B) GAATG. . . . T C GCTG 

Figure 2. Conservation of tlie fHbpRd2R primer binding site. 

Single dots denote bases matching that of the PCR primer (Top 
sequence, 3'-5' direction). Sequence A is representative of 438/443 
(98.8%) meningococcal isolates aligned (Table S2). Sequence B features 
the fHbpRd2R primer site mismatch. This was found in two isolates 
within the MGL (0.22%) and 25% of N. lactamica sequences aligned. 
doi:1 0.1 371/journal.pone.0089921 .g002 

template for the round two reaction appears to have compensated 
for the mismatch as PCR products of appropriate size were 
obtained for all extracts. 

Discussion 

For the purposes of enhanced surveillance, the genotyping of 
meningococcal Factor H-Binding Protein (fHbp) from non-culture 
cUnical specimens is required. To this end, this study describes the 
development and validation of a sensitive PCR-based assay to aid 
vaccine coverage predictions and provide ongoing epidemiological 
surveillance data. 

To measure the sensitivit)' of the PCR assay, analytical 
sensitivity experiments were performed using diluted DNA 
extracts from meningococcal isolates. The use of cultured isolates 
allowed for direct quantification of pure meningococcal DNA. In 
contrast, the use of clinical specimens may have been complicated 
by the presence of human DNA within the extracts. For the nested 
PCR protocol, an analytical sensitivity limit of 6 fg/ |J,L of DNA 
(<3 genome copies/|a.L) was calculated when using 5 |J,L of DNA 
extract in the PCR round-one reaction. At lower DNA concen- 
trations, stochastic DNA sampling effects come into play and the 
likelihood of obtaining an amplified product becomes dependent 
on the volume of extract used in the round one reaction. This is 
illustrated in Figure 1 , in which the number of extracts producing 
ampHcons increases as the round one DNA extract volume used is 
increased. Based on these data, the assay's sensitivity can be 
increased by a factor of 10 (to 600 ag/|iL) by using 20 |xL of 
extract in the first PCR round (in a 100 [ih reaction). Whilst 
maximum sensitivit}' is desired, the assay is to be used routinely 
and the adoption of a 100 |xL round one reaction per specimen 
would be impractical not only in terms of use of consumables but, 
more importantly, in the consumption of the valuable DNA 
extract. A 10 |xL/50 |lL reaction volume would be more 
appropriate and, based on these data, should amplify the majority 
of specimens containing a DNA concentration of ^600 ag/|tL. It 
is recommended that this reaction volume be used in the first 
instance. Increasing the extract volume is of course one of several 
parameters that could be amended in order to gain a product from 
specimens that prove difficult to amplify. These also include 
increasing the PCR cycle number. Although the assay has been 
designed to use as a low cycle number as possible to prevent the 
production of non-specific amplicons, increasing the first round 
cycle number could provide additional sensitivity if required. 
Furthermore, increasing the volume of round one product added 
to the second PCR round reaction has been shown to increase the 
final PCR product concentration. 

The dilution series was tested using the cfr.4-directed TaqMan 
real-time PCR assay used for routine case confirmation. The 
results were comparable to those of the nested-PCR analytical 
sensitivity, indicating that, using this assay, amplification of JHbp 
should be possible for all specimens with meningococcal DNA 



concentrations of S6 fg/^L. This corresponds to an approximate 
Ct value of 38 in this TaqMan assay. Of the non-culture 
specimens received by the MRU in 2011-2012, approximately 
94% achieved a Ct value of ^38. This proportion of clinical 
specimens should therefore yield a PCR product when using at 
least 5 (xL of extract in PCR round one. By increasing the PCR 
round one extract volume to 20 jlL, this predicted proportion 
could be increased to 98% (<Ct41). 

The ability of the assay to amplify ^/7/A/7 from the widest possible 
array of meningococcal specimens is of utmost importance. The 
conservation of the primer binding sequences was assessed by 
testing the assay against a large panel of meningococci, 
representative of major disease-causing clonal populations in the 
UK and other westernised countries. By pairing 96 clinical 
specimens with cultured meningococcal isolates from the same 
disease cases, the results attained could be compared to the data 
produced using whole genome analysis of the isolates to confirm 
the accuracy of the assay and identify any discrepant results. 

PCR products were obtained from 86.5% of the clinical 
specimens on the first attempt. After repeated testing with the 
nested PCR protocol, a total of 91 of the 96 specimen extracts 
yielded products, highlighting the benefit of repeat testing of 
unamphfied extracts and making amendments to the PCR 
parameters, when necessary. For the 96 corresponding meningo- 
coccal isolates, the round one and two PCR protocols were tested 
individually. During PCR round two, one of the isolate extracts 
failed to amplify. This was later found to be due to mismatches in 
the reverse primer binding site. This isolate appears to be one of 
two isolates within the MGL with this identical discordant 
sequence at this site (Figure 2). The two isolates are not of the 
same clonal complex (ST-282 CC and ST-32 GC) and die 
polymorphic site is within the flanking gene (NMB1871), 
suggesting it is simply a rare variation at this site. Interestingly, 
the extract from the clinical specimen matching this validation 
isolate (Mil 909694) was successfully amplified after both nested 
PCR rounds. Subsequent sequencing of this amplicon showed that 
the JHbp allele matches that of the corresponding isolate (allele 69). 
This is a rare allele in the UK and it is associated with the similarly 
rare ST-282 complex, that of the isolate. The JHhp allele match is 
therefore unlikely to have resulted from cross contamination from 
another specimen/extract. In this case, the increase in the 
template DNA concentration during PCR round one may have 
improved the likelihood of semi-complementary primer binding in 
the second round. Early round two products would therefore 
feature the f.xdct JHhpRd2R primer site, and would more efficiently 
amplify from that point. If this is indeed the case, the use of a 
nested PCR procedure may compensate for this rare second round 
primer mismatch and perhaps other primer disparities that may be 
encountered. 

The sequencing results illustrate the conservation of the 
sequencing primers to sites flanking/ within a range of JHhp alleles. 
For the isolates, all amplicons contained JHbp alleles matching 
those assigned from previously acquired whole genome data. All of 
the PCR products from the specimens were sequenced successfully 
and matched their corresponding isolate in terms of JHbp allele. 
InitiaUy, two specimens (Mil 922854 & Mil 915746) produced 
alleles differing from the matched isolate by one nucleotide base, 
however, repeated PCR and sequencing of the extracts resolved 
the discrepancies. The resulting JHbp alleles matched those of the 
respective isolates suggesting PCR copying errors had initially 
occurred. Copying errors were only observed in the clinical 
specimen subset indicating that the low quantity of DNA may be 
increasing the PCR copying error-rate. Copying errors that occur 
during the early PCR cycles may be represented to a significant 
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d<;gT(x- ill the sequence traces of extracts of low starting template 
number. The future adoption of high fidelity 'proof reading' 
polymerase enzymes may be warranted for sequencing of clinical 
specimens. In light of these results, retesting of specimens 
producing JHbp alleles that are not assigned in the PubMLST 
database may also be prudent. 

hAn assessment of primer binding site conservation beyond the 
species level was made applying the assay to isolates of jV. lactamica. 
The closely related commensal does not possess JHbp however 
both sets of PGR primers target neighbouring genes that are 
common to both species, therefore, amplification should occur if 
the sites are suflficientiy conserved. The round one PGR primers 
amplified the locus in all six isolates tested. These findings are 
reflected in a previous study in which the same round one primers 
(1869-2F and 1871Ralt) were used to confirm the absence of JHbp 
in an additional 44 lactamica isolates [26] . The round two PGR 
primer sites were, however, less conserved and the identical 
JHhpRd2R primer mismatch seen in two meningococcal isolates 
(Fig. 2) was shown to be present in 25% of the J^. lactamica genomes 
aligned (Table S2). This finding suggests that the presence of this 
polymorphism in a tiny proportion of meningococcal isolates may 
be a result of recombination with jV. lactamica. Amplification of 
JHbp from all 6 jNI lactamica isolates was achieved when using a 
nested PGR method, providing more evidence to support the 
suggestion that the nested-PGR can compensate for round two 
primer site discrepancies. 

To conclude, this assay has the sensitivity to type JHbp from the 
vast majority of non-culture IMD specimens from England and 
Wales. The development of this assay provides an opportunity to 
complete the epidemiological picture of this important vaccine 
antigen. 

Supporting Information 

Figure SI eBURST and SplitsTree diagrams illustrating the 
distribution of the validation panel isolates among the known 
meningococcal population diversity. eBURST analysis was 
performed using all MLST allelic profile data obtained from the 
PubMLST database on 30/12/2013. Briefly, the eBURST 
program divides STs into groups based on the relatedness of the 
allelic variants at the seven MLST loci. For this analysis, the most 
stringent level was used in which a minimum of 6 identical loci is 
required for group definition. Single locus variants (SLVs, i.e. STs 
that differ at only one of the seven MLST loci) are grouped 
together around a 'founder ST', which possesses the most related 
SLVs. These are predicted to represent an ancestral ST from 
which all other ST variants were derived. Each group roughly 
correlates with established clonal complexes in which constituent 
STs have ^4 loci in common with the ancestral ST. The 
eBURSTs are annotated with the STs of the validation panel 
isolates (highlighted in red) to illustrate their genetic diversity. The 
eBURSTs roughly correlate with 7 predominant clonal complex- 
es: ST-41/44 complex (A), ST-269 complex (B), ST-32 complex 
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